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reviews.  The  report  was  technically  reviewed  by  Dr.  Langway. 
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of  the  use  of  such  commercial  products. 
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THE  TEMPERATURE  STRUCTURE  OF  A  MID-LATITUDE,  DIMICTIC 
LAKE  DURING  FREEZING,  ICE  COVER,  AND  THAWING 

by 

W.H.  Parrott  and  Captain  W.M.  Fleming 


INTRODUCTION 


Vertical  temperature  determinations  in  lakes  are  generally  the  first  and  most  basic  measure¬ 
ments  of  any  limnologic  study.  The  first  series  of  vertical  temperature  profiles  was  made  on 
Scottish  lakes  by  Jardine  during  1812-1814  (Hutchinson,  1957).  Thoreau  was  the  first  American 
to  write  about  thermal  stratification  when  he  lived  at  Walden  Pond  in  the  ISOO’s,  but  he  did  not 
take  enough  readings  to  perceive  the  shape  of  the  thermal  curve  (Deevey,  19-12).  The  work  of  Birge 
(1910,  1915)  from  1895  to  1935  on  Wisconsin  lakes  marked  the  real  beginning  of  limnologic 
studies  (including  temperature  measurements)  in  the  United  States. 

Although  the  temperature  profiles  of  many  of  the  world’s  temperate  lakes  were  generally  de¬ 
scribed  in  the  last  century,  few  investigators  have  made  continuous  and  detailed  measurements  of 
the  thermal  profile,  especially  during  the  critical  times  of  freezing  and  thawing.  The  logistical 
problems  associated  with  temperature  measurements  while  an  init  ial  ice  cover  is  forming  and  while 
it  is  melting  are  important  reasons  for  neglecting  this  area  of  limnological  research.  The  major 
objective  of  this  research  is  to  describe  the  thermal  relationships  at  these  critical  times. 

Continuous  water  temperature  measurements  were  made  in  a  shallow  Like  in  Michigan  (Bilcllo, 
1968)  during  ice  formation,  growth,  and  decay,  but  the  2-m  lake  depth  precluded  the  development  of 
a  classical  three-layered  stratification.  Other  investigators  periodically  (at  weekly  or  longer  inter¬ 
vals)  obtained  water  temperature  profiles  and  calculated  the  annual  heat  budget  for  a  lake  (Birge, 
1915;  Hutchinson,  1937;  Anderson,  1964;  and  many  others). 

Many  analytical  techniques  have  been  proposed  for  calculating  or  estimating  water  temperature 
profiles  for  various  meteorological  conditions.  Some  investigators  concentrated  on  the  exchange 
of  heat  across  the  air /water  interface  (Edinger  el  al.,  1968;  Dutton  and  Bryson,  1902),  while  others 
developed  theories  for  the  vertical  development  of  the  therniocline  over  time  in  a  deep  lake  (Dake 
and  Harleman,  1969).  Important  parameters  of  heat  loss  in  the  overall  budget,  such  as  nocturnal 
radiation  loss  and  its  effect  on  convection  and  thermocline  development,  were  studied  by  Ander¬ 
son  (1968).  However,  few  theoretical  studies  have  been  verified  try  actual  continuous  measure¬ 
ments  of  lake  temperatures. 

During  the  course  of  this  lake  study,  temperatures  were  collected  at  24  data  points  in  the  lake 
every  hour  on  the  hour  over  an  approximately  two-year  period.  This  resulted  in  a  very  large  popu¬ 
lation  of  data.  Obviously  not  all  of  these  data  were  analyzed  or  reported  here.  However,  all  of 
the  data  are  available  in  the  CRREL  library  to  the  CRREL  scientific  staff  or  outside  members  of 
the  scientific  community  for  further  analysis. 
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STUDY  LAKE 

rin*  l:i  ki  *  chosen  t  ( >r  this  .studs'  is  Post  Pond,  a  small  (Hi.  <i  hectares),  mid-lilt  it  twit*  laki*  located 
Is  kin  iioitli  of  Hanover,  New  liampslnrt*.  at  an  fLcvation  of  kill  m  aliovt*  sea  level  (I'm.  1).  It  is 
pan  of  the  C'oiineet n  ut  River  drainage  system.  Post  Pond  lias  a  maximum  depth  of  11.7  in,  is 
relatively  flat  hottomed  (Fig.  2),  and  is  dimietie  (generally  em  ulates  completely  twice  a  year), 
l.imnologieally,  Post  Pond  is  a  Like  rather  than  a  pond*  because  of  its  depth  (average  6.8  til), 
small  littoral  zone,  and  pronounced  three -layered  thermal  stratification. 


Figure  1.  Locution  ( above)  and  watershed  of  Post  Pond. 

Origin  and  structure 

Post  Pond  is  glacial  in  origin  and  occupies  a  basin  that  was  probably  a  large  kettle.  On  two 
sides  of  the  lake  is  an  inegiilar,  glacially  scoured  landscap*  consisting  of  hills  with  a  maximum 
relict  of  :>SS  m.  The  hills  ;ire  com|x>sed  of  early  Paleozoic  metasodiments  (Post  Pond  is  a  volcanic 
member  of  the  Orfordville  formation)  and  are  coveted  thinly  with  glacial  till,  although  bedrock 
outcrops  on  ridge  crests.  On  the  south  is  a  smaller  ridge  of  glacially  scoured  bedrock  and  morainal 
debus,  while  on  the  north  lies  a  swampy  lowland  in  which  both  inlet  (Trout  Brook)  and  outlet  (Clay 
Brook)  are  located.  Clay  Brook  is  incised  In  glacial  clays  and  silts,  indicating  Ihe'fot nier  prosenet 
of  a  huger  glacial  lake. 

*  A  pond  is  defined  bv  Reid  (1961)  as  a  smalt  body  ol  standing  water  shallow  enough  to  have  looted  plaids 
growing  across  it. 
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Figure  2.  Bathymetric  map  of  Post  Pond.  Contour  interval  =  lm. 

Climate  and  hydrology 

The  climate  of  the  upper  Connecticut  Valiev  is  classified  as  humid  continental  with  short,  cool 
summers  and  long,  cold  winters  (Trewartha,  1 954 ).  Sudden  and  rapid  chanties  in  the  weather,  a 
large  daily  and  annual  range  of  temperature,  and  an  equable  yearly  precipitation  distribution  are 
characteristics  of  this  climate.  According  to  the  climatological  records  of  the  Knvironmental  Science 
Services  Administration  -  Weather  Bureau  (1966),  the  average  July  temperature  for  Hanover  is  20.6  C, 
and  the  average  January  temperature  is  -6.9  C.  The  average  yearly  temperature  is  7.1  C.  Annual 
precipitation  is  about  94  cm,  one-quarter  to  one-third  of  which  is  snow. 


The  lake  drains  a  watershed  of  36.1  km2  which  extends  to  the  east  and  north  of  the  lake  (Fig. 
1).  The  major  inlet  stream  'Trout  Brook)  supplies  an  estimated  90' o  of  the  available  runoff  water 
for  the  lake,  and  originates  near  Smarts  Mountain  (9S2-m  elevation)  at  the  eastern  boundaty  of  the 
watershed.  The  outlet  stream  (Clay  Brook)  is  near  the  inlet  stream  on  the  north  side  of  the  Like 
(Fig.  3),  and  flows  northwest  for  5  km  into  the  Connecticut  River.  Discharge  measurements  made 
at  low  ‘low  periods  in  the  fall  of  1968  with  a  Gurley  Current  Meter  showed  that  the  inlet  discharge 
(0.03  mVsec)  was  more  than  twice  the  outlet  discharge.  Flooding  occurs  yearly  in  April  due  to 
spring  snowmelt  (Fig.  4),  and  Trout  Brook  flows  at  an  estimated  10  to  15  times  its  average  dis¬ 
charge  at  that  time. 


Figure  3.  Aerial  view  o(  Post  Pond. 

Local  .mound water  levels  around  the  lake  axe  high  (within  1.5  tn  of  the  ground  surface)  and 
fluctuate  witli  spring  snowmelt  and  heavy  rains.  Shallow  groundwater  aquifers  probably  account 
for  a  significant  amount  of  the  subsurface  inflow  to  the  lake.  Local  groundwater  seepage  into  the 
lake  can  occur  during  the  winter  since  the  snow  covered  soils  remain  essentially  unfrozen  (Likens 
et  ill.,  1967).  The  contribution  of  deep  springs  to  the  hvdrology  of  the  lake  is  unknown. 

An  investigation  of  the  ice  texture,  ice  structure,  and  growth  rate  of  lake  ice  on  Post  Pond  was 
made  during  the  winter  of  1956-57  (Ragle,  1963). 

Organic  productivity 


The  lake  has  a  relatively  high  rate  of  organic  productivity  and  can  be  termed  mildly  eulroplrie 
by  New  Hampshire  standards  according  to  Likens  (personal  communication,  196S).  Many  species 
of  fish  are  present;  among  the  more  abundant  are  txook  trout,  rainbow  trout,  small  mouth  bass,  rock 
bass,  yellow  perch,  and  pickerel. 
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Figure  4.  Spring  flooding  at  Post  Pond. 


METHODS  AND  PROCEDURES 


The  lake  was  instrumented  over  its  maximum  depth  (11.7  m)  with  a  string  of  24  thermocouples 
(Fig.  5).  The  vertical  spacing  of  the  thermocouples  is  presented  in  Figure  7.  The  thermocouples 
were  anchored  in  the  bottom  sediments,  attached  to  a  nylon  line,  and  positioned  at  the  surface  of 
the  lake  by  a  Styrofoam  spar  float.  The  top  six  thermocouples  were  spaced  with  respect  to  the 
surface  of  the  lake,  and  the  lower  eighteen  thermocouples  were  spaced  with  respect  to  the  bottom 
of  the  lake.  One  thermocouple  was  emplaced  25  cm  into  the  bottom  sediments  at  the  end  of  a 
cone-shaped  anchor.  The  lead  wire  extended  along  the  lake  bottom  and  terminated  at  the  instru¬ 
ment  wanigan  located  on  the  shore.  In  addition,  continuous  air  temperatures  and  water  tempera¬ 
tures  of  the  incoming  steam  were  recorded.  Temperature  measurements  were  accurate  to  ±0.1°C. 


Temperature  information  was  collected  from  the  thermocouples  by  an  automatic  data  acquisi¬ 
tion  system  (Model  2010D  Dymec)  with  a  paper  punch  tape  output.  The  major  components  were  a 
digital  voltmeter  and  a  Joseph  Kaye  Reference  Junction.  During  the  first  portion  of  the  season  a 
0°C  temperature  reference  junction  was  used,  but  as  the  lake  approached  freezing  the  signal  dimin¬ 
ished  and  a  150°C  reference  junction  was  substituted.  To  obtain  diurnal  detail,  a  1-hour  sampling 
interval  was  used.  A  computer  program  was  designed  to  convert  the  paper  punch  tape  data  recorded 
in  millivolts  to  temperature  in  degrees  Celsius. 


Other  meteorological  parameters  continuously  recorded  included:  wind  direction  and  velocity 
(Scientific  Associates  Anemometer  No.  470),  incoming  shortwave  solar  radiation  (Kahl  Radiometer 
No.  28  AM  100),  and  net  all-wave  solar  radiation  (Kahl  Radiometer  No.  29  AM  100)  over  the  lake  sur¬ 
face.  Infrequent  observations  and  measurements  included:  precipitation,  relative  humidity,  water 
level  fluctuations  of  the  lake,  snow  depth,  snow  density,  and  ice  thickness. 
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Figure  5.  Instrumentation  schematic  of  Post  Pond. 

lVmpt'iature  measurements  within  the  bottom  sediments  were  measured  with  a  thermistor  probe 
(designed  hv  Likens  and  Johnson,  1969).  The  thermistor  was  read  with  a  modified  Wheatstone 
bridge  which  gave  an  accuracy  within  f0.1°C. 

Water  samples  were  taken  with  a  nonmetallic  water  sampling  bottle  during  various  regimes 
throughout  the  season.  Conductivity  and  cation  analyses  were  made  in  the  CIv!?EL  aqueous  chem¬ 
istry  laboratory. 

Samples  of  winter  ice  were  collected  for  crystallographic  analysis.  A  Whitney  Underwater 
Daylight  Meter  was  used  to  measure  the  transmission  of  solar  radiation  through  the  water. 

RESULTS  AND  DISCUSSION 

Post  Pond  tias  three  pronounced  thermal  layers  at  the  time  of  summer  stratification.  The  upper 
layer  of  more  or  less  uniformly  warm,  circulating,  and  sometimes  turbulent  water  is  termed  the 
epilimmon  (Hutchinson,  1957).  The  lower  layer,  deep,  cold  and  relatively  undisturbed,  is  called 
the  hypolimnion  (Birge,  1910).  In  between  is  a  region  of  rapid  decrease  in  temperature,  commonly 
known  as  the  metalimnion  or  thermocline.  Hutchinson  defined  the  therntocline  as  the  zone  of  maxi¬ 
mum  rate  of  decrease  in  temperature.  Convective  mixing  between  the  epilimnion  and  hypolimnion 
is  essentially  nonexistent.  At  times  of  spring  and  autumn  overturn  (circulation)  the  lake  mixes 
completely.  Therefore,  during  summer  and  autumn  when  the  lake  is  stratified,  the  thermocline  acts 
as  a  barrier  to  mixing  between  upper  (epilimnion)  and  lower  (hypolimnion)  water. 

Summer  stratification 

Hourly  readings  of  lake  temperature  began  on  27  August  1968,  when  the  lake  was  well  strati¬ 
fied  and  in  a  condition  of  extreme  thermal  stability  (Fig.  6).  The  epilimnion  was  isothermal  at 
20.6°C  to  a  depth  of  4.5  m.  As  average  daily  air  temperatures  gradually  decreased,  the  epilimnion 
deepened  to  5.1  m  at  an  isothermal  temperature  of  18.4  C  on  15  September. 

Secondary  and  tertiary  thermocline  development  in  the  epilimnion.  On  16  September  a  warming 
period  began  that  lasted  approximately  15  days  and  resulted  in  the  development  of  secondary  and 
tertiary  t tiermociines  in  the  epilimnion.  The  warming  period  was  started  by  clear  skies  and  little 
wind,  combined  with  increased  air  temperatures. 

Tlie  differences  in  the  temperature  profiles  of  the  epilimnion  between  15  September  and  21 
September  are  striking  (Fig.  7).  On  15  September  the  entire  epilimnion  was  isothermal  at  18.4°C, 
and  the  temperature  curve  was  classical,  with  one  primary  thermocline  extending  from  the  base  of 
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Figure  7.  Secondary  and  tertiary  thermocline  development. 


the  epilimnion  (5.1  m)  to  the  top  of  the  hypolimnion  (8.7  m).  The  primary  thermocline  exhibited  a 
rate  of  temperature  decrease  of  2.7°C/m.  The  primary  thermocline  was  completely  unaffected  by 
subsequent  changes  in  the  epilimnion  in  the  following  6  days. 

Warm  air  temperatures  and  negligible  wind  on  16  September  resulted  in  the  development  of  a 
secondary  thermocline  in  the  epilimnion,  which  extended  from  0.8  m  to  2.7  m  (Fig.  7)  and  showed 
a  rate  of  temperature  decrease  of  0.8°C/m.  The  primary  thermocline  was  identical  to  the  one 
described  for  15  September. 

Continued  high  air  temperatures  combined  with  clear  skies  and  little  wind  resulted  in  the 
development  of  a  tertiary  thermocline  on  21  September  (Fig.  7).  Ttiis  small  thermocline  extended 
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from  0.2  m  to  0.5  hi  and  exhibited  a  rate  ol'  temperature  decrease  of  6.0'  C/ni.  A  secondary  tlier mo¬ 
dule  that  extended  from  0.5  in  to  5.1  in  and  showed  a  rate  of  temperature  decrease  of  0.5°C/m  was 
also  present.  Attain,  the  primary  thermoeline  was  identical  to  the  one  of  the  previous  6  days. 

End  of  the  summer  strut  ideation  period.  The  warming  period  ended  on  1  October,  and  the  sec¬ 
ondary  and  tertiary  thermo, dines  disappeared  as  colder  air  temperatures  and  wind  mixed  the  entire 
epilimnion  isothermally  at  17.6  0  to  a  depth  of  5.7  m  (Fig.  6).  This  date  can  be  termed  the  end  of 
the  summer  stratification  period  and  the  beginning  of  autumnal  mixing. 

Autumnal  mixing  and  thermoeline  disappearance 

Continued  deepening  of  the  primary  thermoeline.  Cool  air  temperatures  and  wind  enabled  the 
thermal  mixing  process  to  cool  the  epilimnion  gradually  and  to  deepen  the  thermoeline  after  1 
October.  By  9  October  the  top  of  the  thermoeline  was  at  6.9  m  with  the  epilimnion  isothermal  at 
14.7  C  (Fig.  6).  On  20  October  the  thermoeline  deepened  to  8.1  m  with  an  isothermal  epilimnion 
at  11.8'C. 

The  thermoeline  was  at  a  depth  of  9.3  m  on  8  November  (with  an  isothermal  epilimnion  of 
9.4  C),  and  by  12  November  at  1000  hours  was  at  9.9  m  (with  an  isothermal  epilimnion  of  8.0°C). 

Autumn  overturn.  The  autumn  overturn  is  defined  as  the  time  during  which  complete  mixing 
occurs  throughout  the  entire  vertical  profile  of  the  lake.  The  overturn  period  began  at  Post  Pond 
on  12  November  at  2200  hours  when  the  entire  lake  became  isothermal  at  7.7°C.  The  water  was 
all  at  the  same  density  and  therefore  extremely  vulnerable  to  mixing. 

As  air  temperatures  decreased,  the  lake  cooled  isothermally  to  7.0  ’C  on  14  November,  to  6.0nC 
on  18  November,  to  5.0nC  on  23  November,  and  to  4.0°C  (temperature  of  maximum  density)  on  29 
November  (Fig.  6).  During  this  time  of  isothermal  cooling  winds  were  light,  but  quite  sufficient 
to  mix  the  lake  completely.  Often  the  top  0.5  m  of  the  lake  cooled  0.5°C  more  than  the  underlying 
11.2  m,  but  within  8  hours  light  winds  were  sufficient  to  mix  the  lake  and  cool  the  entire  12-m 
column. 

TEMPERATURE,  "C 


Figure  8.  Winter  temperature  profiles. 
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Reverse  stratification  was  evident  on  1  December,  with  the  upper  0.5  m  colder  than  3  C  and 
the  bottom  remaining  at  3.7°C  (Fig.  8).  By  6  December  light  winds  were  sufficient  to  circulate 
the  lake  for  the  last  time  and  mix  it  completely  to  3.2°C.  Calm  conditions  and  freezing  air  tempera¬ 
tures  led  to  the  development  of  permanent  reverse  stratification  on  7  December,  signifying  the  end 
of  the  autumn  overturn.  The  overturn  had  lasted  from  12  November  to  7  December  -  a  total  of  25 
days.  An  important  cause  of  the  long  overturn  period  was  the  absence  of  storms  or  strong  winds, 
which,  if  present,  could  have  combined  with  cold  air  temperatur  es  to  bring  the  lake  to  4  C  in  a 
few  days. 

Winter  period  of  ice  cover 

Establishment  ol  ice  cover.  On  the  night  of  8  December  calm,  cold  conditions  led  to  the  forma¬ 
tion  of  an  initial  ice  skim  which  melted  the  following  day.  Similar  conditions  the  following  night 
led  to  the  initiation  of  a  permanent  ice  cover.  Rapid  growth  of  clear  ice  resulted  in  a  cover  approxi¬ 
mately  4  cm  thick  the  following  day.  The  resulting  inverse  stratification  on  10  December  intensi¬ 
fied  (Fig.  8). 

Ice  covered  period  -  cooling  ol  the  water  mass.  The  ice  cover  lasted  from  10  December  1968 
to  19  April  1969  (130  days,  or  36%  of  a  year).  During  this  time  the  ice  cover  thickened  to  a  maxi¬ 
mum  of  62  cm  on  7  March  (42  cm  of  clear  ice,  20  cm  of  snow-ice*),  but  the  temperature  structure 
of  the  water  changed  only  slightly.  During  March  1968,  ice  samples  were  taken  from  six  sites  on 
the  pond.  Some  of  these  data  are  presented  in  Appendix  A.  The  difference  between  the  snow-ice 
and  the  clear  ice  is  quite  discernible.  Profiles  from  the  four  months  of  ice  cover  show  that  the 
upper  7  m  cooled  slightly  while  the  lower  5  m  warmed  slightly  (Fig.  8). 

Ice  covered  period  -  heat  from  bottom  sediments.  Slow  warming  of  the  lower  5  nt  of  water  was 
evident  (see  Fig.  8).  By  15  February  the  lower  part  of  the  hypolimnion  had  warmed  to  3.5  C  (from 
3.2°C  on  10  December);  on  5  March  it  had  warmed  to  3.8°C;  and  by  15  March  it  had  warmed  to 
4.0°C.  This  heat  was  furnished  to  the  lake  from  bottom  sediments  in  the  form  of  solar  heat  stored 
during  the  summer  months,  to  be  released  to  the  water  during  the  period  of  ice  cover.  Although  a 
small  amount  of  visible  solar  radiation  may  penetrate  2  m  of  lake  ice  and  affect  ice  growth  at  the 
ice/water  interface  (Swinzow,  1966),  no  infrared  radiation  can  be  expected  to  penetrate. 

Results  with  a  Whitney  Underwater  Daylight  Meter  on  8  January  (with  a  38-cnt  ice  cover)  showed 
that  approximately  10%  of  total  surface  light  penetrated  1  nt  of  ice  and  water,  while  only  about 
0.1%  penetrated  5  m  (see  Fig.  9).  In  addition,  the  ice  was  continually  shielded  from  solar  radia¬ 
tion  by  a  substantial  snow  cover  (average  depth  25  cm)  from  22  December  to  18  March.  This  reduced 
the  penetration  more. 

Thermal  gradients  in  bottom  sediments  were  significant  portions  of  the  total  heat  budget  of 
lakes  according  to  Johnson  and  Likens  (1967)  and  Likens  and  Ragotzkie  (1965).  Studies  of  Likens 
and  Johnson  (1969)  showed  that  the  bottom  sediments  of  two  small  Wisconsin  lakes  contributed 
10-12%  of  the  total  annual  heat  budgets. 

Thermal  gradients  were  measured  on  13  March  (approximately  the  time  of  maximum  ice  thick¬ 
ness)  to  a  depth  of  9  m  in  the  sediments  (Fig.  10).  An  average  thermal  gradient  of  about  0.5  C/m 
was  found  for  the  upper  3  m  of  sediments  and  a  steady-state  thermal  gradient  of  about  0.07  C/nt 

*  Snow-ice  is  a  frozen  slush  layer  that  forms  over  the  clear  ice  from  a  mixture  ol  snow  and  lake  water.  The 
accumulating  load  of  fallen  snow  causes  the  clear  ice  to  fracture;  water  ascends  through  the  fracture 
pattern,  mixes  with  the  snow,  and  freezes.  The  contact  between  clear  ice  and  snow-ice  is  easily  dis¬ 
cernible  (see  App  A). 
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Figure  9.  Depth  of  ice  and  water  v's  percentage  of  light  penetration. 
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Figure  10.  Temperature  profile  on  13  March  1969. 
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was  found  for  the  underlying  6  m.  The  steady-state  thermal  gradient  is  in  part  a  measure  of  the 
terrestrial  geothermal  gradient,  and  other  studies  have  calculated  that  the  geothermal  #adient 
accounts  for  about  half  of  the  steady-state  gradient  measured  in  lake  sediments  (Johnson  and 
Likens,  1967).  A  steady-state  thermal  gradient  of  0.20°C/m  was  measured  in  the  bottom  sediments 
of  a  small  Wisconsin  lake  (Likens  and  Johnson,  1969);  this  is  approximately  three  times  the  grad¬ 
ient  measured  for  the  sediments  of  Post  Pond. 

Considering  only  the  larger  thermal  gradient  of  0.5°C/m  in  the  upper  3  m  of  sediments,  the 
heat  flux  through  the  sediments  to  the  water  can  be  calculated  according  to  the  general  heat  con¬ 
duction  equation: 

T 

<t>  =  k  A- 
L 

where: 

=  heat  flux  through  sediments  (cal/cm2  day) 

=  thermal  conductivity  of  sediments  (cal  cm/cm2  sec°C) 

=  average  thermal  gradient  in  sediments  (°C/cm). 

Since  the  bottom  sediments  of  Post  Pond  are  similar  to  those  of  Tub  Lake,  Wisconsin,  a  simi¬ 
lar  thermal  conductivity  (1.36  x  1(TS  cal  cm/cm2  sec°C)  is  assumed  (Likens  and  Johnson,  1969; 
McGaw,  in  prep).  Thus,  the  heat  flux  conducted  through  the  upper  3  m  of  sediments  into  Post 
Pond  during  maximum  ice  cover  was  -6.8  x  10“‘  cal/cm2  sec  or  -0.59  cal/cm2  day. 

The  amount  of  heat  conducted  into  the  lower  5  m  of  water  is  enough  to  account  for  the  gradual 
increase  in  heat  content.  The  amount  of  heat  introduced  during  the  period  10  December  to  15  March 
is  calculated  according  to; 

D, 

Qw  =  S  (Ttl-  Ttg)  (AJ»  (fw)  (pcp) 

—i 

where: 

Qw  =  heat  content  of  the  water,  cal/cm3 
Dj  =  depth  of  6.7  m 
Dg  =  depth  of  11.7  m 
Ttl  =  water  temperature  on  10  Dec,°C 

T t2  =  water  temperature  on  15  Mar,  °C 

A/i  =  thickness  of  water  layer  at  each  depth  interval,  cm 
f  =  weighting  factor  from  hypsographic  curve  for  relative  area  at  each  depth 
cp  =  specific  heat  capacity  of  water,  cal/g"C 
p  =  density  of  water,  g/cms. 

Thus,  the  amount  of  heat  gained  by  the  lower  5  m  of  water  during  this  period  was  approximately 
51.5  cal/cm3. 

The  amount  of  heat  conducted  from  the  bottom  sediments  during  the  period  10  December  to 
15  March  is  calculated  according  to: 

<?R  = 


<t>  \t 
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Qs  heat  flux  m  upper  3  m  of  sediments, cal/cinJ 
<l>  heat  flux  through  sediments,  cal/em1  day 
\t  ice  covered  period,  in  days. 

Tlie  heat  flux  in  the  sediments  near  the  end  of  t he  iee-cwered  period  (-0.59  eal/cm'  day  on 
13  Marcti)  is  used,  as  it  represents  a  smaller  gradient  than  would  be  expected  at  the  beginning  of 
the  winter  (Likens  and  Johnson,  1969),  and  thus  is  a  conservative  measure  of  the  heat  available 
from  the  sediments: 

Qs  (-0.59)  (96)  -56.6  eal/eniL 

The  nearly  equal  values  of  Qw  and  Qs  support  the  probability  that  heat  from  bottom  sediments  is 
enough  to  account  for  warming  the  water  mass  during  ice  cover.  These  results  support  Bilello's 
(1968)  conclusion  that  bottom  muds  constitute  the  major  heat  source  during  ice  cover.  These 
results  suggest  that  little  heat  was  udvected  into  Post  Pond  by  deep  groundwater  sources  during 
the  winter. 

These  results  must  be  interpreted  as  having  validity  only  for  the  sediments  in  the  deepest  part 
of  the  lake,  as  Likens  and  Johnson  (1969)  show  that  sediments  near  the  shoreline  of  lakes  are 
thermally  much  more  active. 

Takagi’s  (1970)  use  of  spline  functions  to  analyze  the  temperatures  observed  in  a  lake  enabled 
us  to  detect  the  period  during  which  only  vertical  heat  conduction  was  the  mechanism  of  heat  trans¬ 
fer.  If  the  scattering  due  to  observation  errors  is  disregarded,  a  quantity  called  ‘the  integral 
residual”  becomes  a  minimum  in  this  period.  The  standard  deviation  of  the  temperature  observa¬ 
tion  can  he  found  from  the  scattering.  Only  the  differential  equation  is  needed  in  this  computa¬ 
tion;  none  of  the  boundary  and  initial  conditions  are  required. 

Ice  covered  period  -  lute  winter  cooling  of  bottom  water.  Alter  warning  gradually  to  4.0’C, 
the  bottom  2  in  of  water  cooled  suddenly  (within  24  hours)  to  3.1°C  o  21  March  (Fig.  11).  Although 
it  occurred  at  all  depths,  a  similar  spring  cooling  period  was  reported  by  Bilello  ( 1968).  The 
bottom  ""iter  cooled  further  to  2.9°C  on  24  March,  and  remained  at  that  temperature  until  27  March 
when  gan  to  warm.  By  29  March  the  water  was  at  3.4°C  and  warmed  gradually  to  3.7  C  on  9 
April,  ihe  short  cooling  period  lasted  about  6  days  (21  March  to  27  March).  Two  possible  explan¬ 
ations  are  offered:  1)  the  bottom  sediments  became  thermally  drained,  thus  causing  the  bottom 
water  to  cool  and  convectively  mix  with  the  upper,  cooler  water;  and  2)  warming  air  temperatures 
(daily  average  air  temperatures  increased  from  1.9°C  to  3.0°C  during  this  period)  caused  snowmelt 
at  higher  elevations  in  the  watershed  to  activate  deep  groundwater  influx  to  lower  levels  of  the 
lake. 

The  first  explanation  is  rejected  because  of  the  rapid  increase  in  water  temperature  after  the 
6-day  cooling  period.  It  is  doubtful  that  the  sediments  would  become  thermally  reactivated  quickly 
after  being  drained.  The  second  explanation  is  more  plausible,  as  snowmelt  activation  of  ground- 
water  aquifers  can  be  a  rapid  consequence  of  increased  air  temperatures  (Fleming,  1966).  Since 
tlie  forest  soils  of  the  lower  elevations  of  New  Hampshire  remain  unfrozen  during  the  winter  (Likens 
et  a/.,  1967),  snowmelt  water  quickly  percolates  to  the  groundwater  table,  which  is  relatively 
high  in  ttie  Post  Pond  watershed.  Further  geohydrological  study  of  the  watershed  is  needed  to 
support  this  interpretation. 
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Figure  11.  Late  winter  cooling  of  bottom  water. 

Melting  of  the  ice  cover  -  melting  process.  The  melting  period  lasted  from  7  March  to  19  April 
(43  days).  Ice  generally  melts  from  the  top  due  to  the  development  of  a  high  positive  net  radiation 
in  the  spring  (Hutchinson,  1957).  Temperature  measurements  at  the  ioe/water  interface  combined 
with  ice  thickness  measurements  at  Post  Pond  indicate  that  the  process  may  be  more  complex, 
especially  with  the  added  parameter  of  snowmelt  runoff. 

The  melt  season  at  Post  Pond  can  be  conveniently  divided  into  three  periods:  1)  initial  radia¬ 
tion  melting  of  the  upper  snow-ice  from  8  March  to  2  April;  2)  melting  of  the  clear  ice  from  under 
the  ice  cover  (2  April  to  10  April);  and  3)  melting  of  the  remainder  of  the  snow- ice  from  10  April 
to  19  April.  Initial  melting  of  the  snow-ice  by  radiation  began  on  8  March  and  continued  slowly 
until  2  April,  melting  the  upper  5  cm  of  the  snow-ice. 

Snowmelt  water  began  to  flow  into  the  lake  via  TYout  Brook  and  several  small  ephemeral 
streams  on  26  March.  Flooding  of  the  lowlands  occurred  from  2  April  to  10  April,  and  it  was  during 
this  peak  snowmelt  period  (8  days)  that  the  entire  clear  ice  section  of  the  cover  (42  cm)  melted 
(Fig.  12  shows  the  growth  and  melt  of  clear  ice).  The  significant  aspect  of  this  melting  phenomenon 
is  that  temperature  and  ice  thickness  measurements  show  that  it  occurred  on  the  underside  of  the 
ice  and  was  mainly  due  to  heat  supplied  by  snowmelt  water.  Temperatures  of  the  influent  stream 
(Trout  Brook)  averaged  3.3°C  during  this  8-day  period.  Radiational  heating  through  15  cm  of  snow- 
ice  was  considered  small  (Lyons  and  Stoiber,  1959;  Thomas,  1963). 

The  lake  began  to  melt  first  around  the  shoreline,  creating  a  moat  around  the  lake  that  enlarged 
in  width  as  melt  continued  (Fig.  13).  Snowmelt  water,  mainly  from  the  east  and  north  sides  of  the 
lake,  contributed  to  the  enlargement  of  the  moat  and  supplied  heat  for  melting  the  clear  ice  layer 
at  the  base  of  the  ice  cover.  The  development  of  drainage  holes,  thought  by  some  investigators 
(Williams,  1968)  to  be  a  major  cause  of  spring  melt,  were  not  observed  on  Post  Pond. 
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Figure  12.  Growth  and  melt  of  clear  ice 
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Figure  13.  Initial  melting  and  moat  formation  around  shoreline  in  spring. 


The  third  and  final  period  of  melt  occurred  from  10  April  to  19  April  and  resulted  in  the  com¬ 
plete  disappearance  of  ice  from  Post  Pond.  Heat  for  melting  the  remaining  15  cm  of  snow-ice  was 
supplied  by  snowmelt  water,  radiation,  and  wind  mixing.  Candling  (melting  at  crystal  boundary) 
of  the  thin  cover  observed  during  this  period  (Fig.  H)  contributed  to  the  disintegration  of  the  cover 

Melting  of  the  ice  cover  -  thermal  instability  of  the  water  mass  during  melt  season.  Thermal 
instability*  of  the  water  mass  was  observed  for  9  days  (11-19  April)  during  the  last  period  of  the 
melt  season  (Fig.  15).  (Thermocouple  readings  were  checked  against  independent  thermistor 
readings  and  found  to  be  accurate.)  Thermal  instability  was  reported  by  oilier  investigators 
(Ruttner,  1952;  Bilello,  1968);  Ruttner  believed  it  to  be  a  result  of  chemical  stratification. 

*  Tlirnnat  instability  is  defined  as  a  |>eriod  during  which  water  at  either  greater  or  less  than  4  C  underlies 
water  at  1  C. 
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Figure  14.  Candling  of  thin  ice  cover  during  spring  melting. 
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Figure  15.  Thermal  instability  and  conductivity  during  melt  season. 


Beginning  on  11  April  the  upper  2.5  m  ol'  water  was  slightly  thermally  unstable.  By  15  April  trie  upper 
meter  of  water  had  warmed  beyond  4°C  (temperature  of  maximum  density),  and  water  from  1.2  m  to 
the  bottom  was  thermally  unstable.  By  19  April  the  unstable  condition  ceased  and  the  entire  lake 
was  isothermal  at  3.5°C.  By  25  April  slight  stratification  began  to  develop,  witli  the  upper  9.5 
m  isothermal  at  4.9 'C. 

During  the  period  of  instability,  thermal  profiles  showed  significant  changes  within  a  few 
hours.  Thermal  profiles  at  1100  and  1400  on  18  April  are  shown  in  Figure  16.  Continuous  ?l„  mal 
adjustment  of  the  upper  8  m  is  evident.  Thermal  instability  occurred  from  l.S  m  to  the  bottom  a 
1100  hours,  but  subsequent  mixing  resulted  in  a  layer  of  unstable  water  from  7.5  to  11.0  in  at 

1400  hours. 
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Figure  16.  Effect  of  temperature  and  dissolved  solids 
on  density,  18  April. 

Conductivity  analyses.  The  conductivities  of  water  samples  collected  at  1-m  vertical 
intervals  are  shown  in  Figure  15.  The  water  samples  were  collected  with  a  nonmetallic  sample 
bottle  and  the  conductivities  were  measured  with  a  Leeds  and  Northrup  conductivity  meter.  The 
change  in  conductivities  from  a  time  of  thermal  instability  (18  April)  to  a  time  of  thermal  stability 
(25  April)  is  marked,  especially  in  the  lower  3  m.  On  18  April  the  upper  2  m  of  water  had  conduc- 
t  i v it  ics  averaging  10  /(iiilios/cm  less  than  the  underlying  6  m,  and  t he  bottom  3  m  of  water  had 
conductivities  a  veracity  08  /unhos/om  more  than  the  overlying  7  m,  indicating  chemical  stratifi¬ 
cation.  Between  19  and  25  April,  wher  the  lake  mixed  thoroughly,  a  significant  change  occurred, 
and  values  were  essentially  the  same  throughout  the  depth  of  the  lake.  All  evidence  of  instability 
disappeared  and  the  lake  was  completely  mixed  to  10  m. 

Dissolved  solids  analyses.  Dissolved  solids  of  water  samples  collected  at  1-m  vertical 
intervals  at  1400  hours  on  18  April  are  shown  in  Figure  16.  Analyses  included  the  major  cations 
(Ca,  Mg.  K,  and  Na);  and  for  this  interpretation  the  sulfate  (SOJ  anion,  although  not  measured  on 
these  samples,  was  assumed  to  be  the  major  anion.  Cation  analyses  are  shown  in  Table  I.  This 
assumption  is  reasonable  considering  the  high  sulfate  content  of  New  Hampshire  waters  (Love, 
i960).  The  dissolved  solids  were  stratified  and  increased  with  depth,  showing  the  greatest  in¬ 
crease  (28  mg/liter)  between  9  m  and  10  m. 

Combined  effect  of  temperature  and  dissolved  solids  on  density.  The  density  of  the  water 
mass  calculated  from  the  temperature  data  at  1400  hours  is  shown  in  Figure  16.  Overall  instability 
occurred  from  7.5  m  to  9.2  m,  and  can  be  termed  thermo-chemical  instability  because  it  combines 
the  effects  of  temperature  and  dissolved  solids.  Water  that  is  apparently  8  mg/liter  lighter  than 
the  water  above  occurs  in  the  unstable  layer.  It  is  concluded  that  while  an  increase  with  depth 
in  dissolved  solids  aids  in  explaining  thermal  density  instability,  it  is  obviously  not  the  whole 
explanation. 

The  period  of  instability  correlates  well  with  the  time  of  peak  snowmelt  runoff.  Particulate 
matter  (suspended  solids)  deposited  try  snowmelt  waters  (not  measured  in  this  analysis)  may  have 
contributed  to  further  instability  especially  at  lower  levels  in  the  lake. 
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Table  1.  Cation  analyses,  18  April  1969  (mg/literk 


Depth  ( m ) 

_  K * 

IVaf 

Ca^' 

Surface 

0.20 

2.30 

0.74 

6. 10 

1 

0.30 

2.00 

0.77 

6.40 

2 

0.30 

3.00 

0.80 

7.00 

3 

0.30 

3. 10 

0.80 

7.00 

4 

0.30 

3.00 

0.88 

7.  10 

0 

0.30 

3.40 

0.90 

7. 10 

6 

0.30 

3.40 

0.89 

7.20 

7 

0.30 

3.80 

0.93 

7.60 

8 

0.30 

4.10 

0.93 

7.00 

9 

0.30 

4.20 

0.93 

7.00 

10 

1.00 

7.60 

1.18 

11.80 

11 

1. 10 

9.90 

1.42 

9.00 

Spring  circulation 

Spring  overturn.  On  19  April  three  significant  events  occurred:  1)  t he  ice  cover  completely 
melted;  2)  the  thermally  unstable  profiles  of  the  previous  9  days  disappeared;  and  3)  the 
spring  overturn  period  began.  The  beginning  of  the  overturn  occurred  when  the  entire  lake  mixed 
completely  and  became  isothermal  at  3.4°C  (Fig.  17).  The  overturn  period  lasted  4  days,  until 
23  April,  when  the  lake  warmed  isothermally  to  4.5°C  and  then  began  to  stratify. 

TEMPERATURE,  “C 


Figure  17.  Spring  circulation  and  thermocline  development. 
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Thennochne  development,  The  mixing  process  in  the  spring  development  of  tin*  thermocline 
was  quit t>  different  from  the  mixing  prix-oss  in  tin*  uut tiiiiu  when  the  tlicrmncliiic  was  gradually  destroyed. 
As  i lie  upper  water  roolfd  in  the  autumn  and  became  niort*  dense,  the  wind  helped  to  mix  it  convec- 
tively  with  lighter  water  beneath.  In  this  ease  wind  and  convection  worked  together  to  mix  the 
lake.  In  the  spring,  however,  the  convection  process  and  wind  were  working  against  each  ot her. 

As  the  upper  water  warmed  due  to  increased  air  temperatures  and  radiation,  it  became  less  dense 
and  tended  to  stay  on  top.  Therefore,  the  only  energy  available  to  mix  the  lake  in  the  spring  was 
wind. 

Spring  winds  mixed  the  lake  to  a  depth  of  8  in  on  2b  April  and  warmed  the  upper  water  to  5.7  'C 
(Fig.  17).  Density  differences  had  become  great  enough  on  27  April  so  that  light  winds  were  able 
to  mix  the  lake  to  only  5  m.  Due  to  a  relatively  windless  period  of  5  days,  mixing  occurred  to  only 
2  m  by  May.  Strong  winds  then  overcame  density  differences  to  mix  the  lake  to  a  depth  of  4.5  m 
on  11  May,  with  classical  three-lavered  stratification  quite  evident  (Fig.  17).  At  this  time,  density 
differences  had  become  great  enough  so  that  even  strong  winds  were  not  sufficient  to  develop  the 
thermoeliiie  further.  The  period  of  spring  circulation  was  essentially  over,  and  the  summer  strati¬ 
fication  period  had  begun. 


SUMMARY  AND  CONCLUSIONS 

The  data  show  that  Post  Pond  is  a  typical  diabetic,  mid-latitude  lake;  based  on  the  data, 
the  following  statements  can  be  made. 

1.  Secondary  and  tertiary  thermocline  development  in  the  epilimnion  occurred  during  short 
warming  periods  in  the  early  autumn.  The  rate  of  temperature  decrease  in  the  tertiary  thermocline 
(b.l)  C  /m)  was  more  than  twice  the  rate  of  decrease  in  the  primary  thermocline. 

2.  The  autumn  overturn  Listed  25  days,  whereas  the  spring  overturn  lasted  only  4  days. 

3.  The  entire  lake  tnixtxl  isothermally  in  the  autumn  to  3.2°C,  almost  one  degree  lower  than 
the  temperature  of  maximum  density. 

4.  Temporal  lire  profiles  during  the  4  months  of  ice  cover  show  that  the  upper  7  m  cooled 
slightly  while  the  lower  5  m  warmed  slightly.  The  lower  5  m  of  water  gained  approximately  51.5 
cal/ ciiP  from  10  December  to  15  March. 

5.  The  heat  flux  conducted  through  the  upper  3  m  of  sediments  into  Post  Pond  during  maximum 
ice  cover  was  approximately  0.6  cal/cm2  day,  or  56.6  cal/cm2  from  10  December  to  15  March.  Heat 
from  the  bottom  sediments  is  enough  to  account  for  the  warming  of  the  water  mass  during  the 
existence  of  a  stable  ice  cover. 

6.  A  steady-state  thermal  gradient  of  0.07  C/m  was  found  for  the  deeper  6  m  of  the  sediments 
underlying  the  lake  during  ice  cover.  This  is  approximately  three  times  the  normal  geothermal 
gradient  of  0.02  'm  quoted  by  Van  Orstrand  ( 1039).  This  is  also  approximately  half  as  much  as 
tlie  thermal  gradient  of  0.20  C/m  measured  in  a  Wisconsin  lake  (Likens  and  Johnson,  1969).  The 
thermal  conductivity  of  the  sediments  was  not  measured. 

7.  Late  winter  cooling  of  bottom  water  under  the  ice  cover  may  be  the  result  of  snowmelt  in 
the  areas  adjacent  to  tiu*  lake  causing  activation  of  groundwater  influx  to  lower  levels  of  the  lake. 

S.  The  clear  ice  portion  of  the  ice  cover  melted  in  S  days;  this  was  caused  primarily  by  heat 
supplied  to  the  lake  by  snowmelt  water.  Melting  occurred  on  the  underside  of  the  ice  sheet. 
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9.  Thermal  instability  of  the  water  mass  persisted  for  9  days  during  peak  snowmelt  runoff; 
this  can  be  partially  explained  by  an  increase  in  dissolved  solids  with  depth. 
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Figure  A3.  Side  view  of  sample  collected  Figure  A  t.  Side  vii'w  of  u  sample  collected 

from  site  2  with  approximately  35  cm  of  ice.  from  site  5  with  over  40  cm  of  ice.  20  cm  of 

15  cm  clear  ice  and  20  cm  snow-ice.  The  which  is  clear  ice.  The  difference  between 

dark  band  is  the  original  snow-ice  and  was  this  sample  and  that  for  site  2  represents  the 
probably  frozen  rapidly.  inhomogeneity  of  ice  on  a  particular  body  of 

water. 
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The  temperature  structure  of  P  st  Pond,  a  small  (46.6  hectares),  mid-latitude,  dimictic 
lake  In  west-central  New  Hampshire,  was  studied  during  autumn, winter  and  spring  of 
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underlying  the  lake  were  measured  with  a  thermistor  probe.  Secondary  and  tertiary 
thermooline  development  in  the  epilimnion  occurred  during  short  wurming  periods  in  the 
early  autumn.  The  autumn  overturn  lasted  25  days,  whereas  the  spring  overturn  lasted 
oniy  4  days.  The  entire  luke  mixed  isothermally  in  the  autumn  to  3 .2*C.  During  the 
period  of  ice  cover,  the  lower  5  m  of  water  gained  approximately  51*5  cal/cm .  ,  which 
was  out  plied  by  stored  heat  in  the  bottom  sediments.  A  steady-state  thermal  gradient  o: 

.  07°c’ra  was  found  for  the  deeper  sediments  underlying  the  lake  during  ice  cover.  Late 
winter  cooling  of  bottom  water  under  the  ice  cover  may  be  the  result  of  snowmelt  in  arei 
adjacent  to  the  lake  causing  activation  of  groundwater  influx.  Melting  of  the  clear 
ice  portion  of  the  ice  cover  was  primarily  the  result  of  heat  supplied  to  the  lake  from 
snowmelt  water,  and  occurred  on  the  underside  of  the  ire  sheet.  Thermal  Instability 
of  the  water  mass  persisted  for  9  days  during  peak  snowmelt  runoff;  this  can  be 
partially  explained  by  an  increase  in  dissolved  solids  with  depth. 
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